^  ■ 


U.D,C.  No.  531.52  :  621.372.4  :  621.391.825  ;  551.501.75 


ROYAL  AIRCRAFT  ESTABLISHMENT 
Technical  Report  N0.64049 
Novemher  I964 


TERMINAL  VELOCITIES  OFtTODOW  DIPO^^  USED  IN 

HIGH  ALTITUDE  /flMlffiSSIJEEi&S  ^  V' 

. //  ■■ 

T,  W.  G,  Ddwson,  M.A.,  A.  Inst.P.  '''(P' 


■TECEhIG.AL  library 

EI.S0  SjC 

ABEEDEEH  PEOViLA  GLGU'.; WO. 
STUiP-Th 


SUMMARY 


Various  theoretical  formulae  for  the  rate-of-fall  of  window  dipoles  at 
altitudes  up  to  100  Km  are  shown  to  be  in  close  agreement  with  the  best  avail¬ 
able  experimental  data.  At  low  altitudes  velocity  is  determined  mainly  by 
conventional  aerodynamic  drag.  At  approximately  20  Km,  Lamb’s  theoxy  of 
viscous  flow  is  required  whilst  at  50-70  Km,  corrections  must  be  applied  for 
’slip-flow’.  Finally,  above  about  75  Km,  ’free  molecule  flow’  theory  is 
needed  to  get  good  agreement  with  experiment. 

For  the  most  accurate  measurements  of  wind  structure  it  is  desirable  to 
employ  dipoles  with  the  smallest  possible  terminal  velocity.  This  means  that 
the  dipoles  should  have  low  density  and  the  smallest  possible  thickness. 
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1  INTRODUCTION 

1  -3 

One  well-established  technique  for  the  measurement  of  winds  at  high 
altitudes  above  the  surface  of  the  earth  is  to  release  into  the  atmosphere,  by 
means  of  a  rocket,  a  cloud  of  small  strips  of  metal  known  either  as  ‘window 
dipoles'  or  'chaff.  These  fall  towards  the  earth  under  the  action  of  gravity 
and  are  also  blown  to  and  fro  by  the  wind.  The  dipoles  are  easily  visible  to 
a  radar  equipment  and  thus,  provided  the  rate-of-fall  is  not  too  great,  the 
magnitude  and  direction  of  high-altitvide  air-streams  can  be  computed  from  range 
and  angle  data  obtained  by  the  radar.  In  a  typical  experiment  a  single  cloud 
might  be  formed  from  one  million  dipoles  with  a  total  weight  of  one  pound. 

This  might  be  released  at  85  Km  altitude  and  would  act  as  a  wind-marker  for  a 
time  of  the  order  of  30  minutes  before  becoming  too  diffuse  to  be  tracked  by 
radar. 

In  a  previous  report^,  a  theory  v/as  developed  for  the  rate-of-fall  of 
dipoles  released  at  various  altitudes.  At  that  time,  however,  the  experimental 
data  was  rather  inaccurate  and  it  was  only  possible  to  check  that  the  theory  and 
experiment  were  in  general  agreement.  Since  then,  the  following  developments 
have  occurred; 

(i)  much  more  precise  measurements  have  been  made^  of  the  rate-of-fall 
at  various  altitudes, 

(ii)  more  accurate  measurements  have  been  made  of  the  density,  pressure 
and  temperature  of  the  upper  atmosphere®, 

(iii)  another  report^,  containing  theoretical  formulae  for  rate-of-fall, 
has  become  available. 

The  purpose,  therefore,  of  this  present  report  is  to  review  the  theory 
(Section  2)  and  then  to  make  an  accurate  comparison  between  theory  and  experi¬ 
ment  in  order  to  discover  any  discrepancies  and  to  determine  how  to  make  the 
dipoles  fall  as  slowly  as  possible. 

2  THEORETICAL  FORMUIAE 

2.1  General  theory 

The  general  equation  of  motion  for  a  dipole  falling  through  still  air  is 


=  mg  -  pV^  A  C 


m 


M. 

dt 


(1) 


If 


In  this  equation  the  LHS  represents  mass  times  acceleration  whilst  the  RHS  is 
the  gravitational  force  less  aerodynaraio  drag.  Boyancy  forces  are  assumed 
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zero  .  Also: 

V  =  downward  velocity 

p  =  air  density  ’ 

A  =  dipole  projected  area 

m  =  dipole  mass 

Cjj  =  the  drag  coefficient  (defined  hy  equation  (l)). 

To  find  the  dipole  terminal  velocity,  set  dV/dt  =  0: 

Tr2  2mg 

=  pACj^ 


(2) 


Thus  to  obtain  V  it 
In  general"'^, 


is  necessary  to  calculate  the  drag  coefficient 

Cjj  is  a  function  of  Reynolds  nmber  R  and  Mach  number  M: 


(3) 

(4) 


where  w  =  dipole  width 

(i  s  viscosity  of  air 
c  =  velocity  of  sound. 

The  Reynolds  number  may  be  written  in  terms  of  any  characteristic  dimension  of 

the  body.  In  this  report,  however,  width  is  chosen  because  it  appears  from  the 

results  given  below  that  most  dipoles  tend  to  fall  with  their  length  in  the 

horizontal  plane  and  in  this  situation  the  various  theories  show  that  the  width, 

rather  than  the  thickness  or  length,  is  the  most  important  dimension  affecting 
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Cp  and  V.  The  velocity  of  sound  at  any  altitude  may  be  computed  from 

0^  =  ^  (5) 

where  y  =  ratio  of  specific  heats  of  air  molecules  =  1.41 
P  =  air  pressure. 
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2,2  Low  altitudes 

At  low  altitudes,  calculation  of  is  a  straightforward  matter;  the 
falling  dipole  causes  a  disturbance  in  the  surrounding  air  which  acts  as  a 
continuous  fluid.  Near  ground  level  it  is  found  that 

M  «  VR 

and 

R  >  1  (6) 

and  Cjj  may  be  found  from  standard  text-books  on  low-speed  aerodynamics. 

If  the  main  surface  of  the  dipole  is  in  the  horizontal  plane  then  experi¬ 
mental  results  can  be  approximated  by  the  relations^: 

Cp  ~  2  at  R  >  20,  (7) 

Cjp  ~  a,t  20  >  R  >  2  ,  (8) 

Taking  A  =  Lw  where  L  is  the  dipole  length  and  substituting  for  in  equation 
(2),  it  is  then  found 

=  (f)  •  (5) 

2 

=  (f)  ■  20  >  R  >  2  .  Oo) 

If  a  flat-plate  has  an  edge  at  right-angles  to  the  flow,  (e,g,  the  long 
edge),  then  the  drag  components  are  skin  friction  and  forebody  drag.  The  latter 
contribution  is,  however,  negligible  for  dipoles  at  low  altitudes.  The  standard 
formula  for  skin  friction  in  laminar  flow  conditions  is  that  due  to  Blasius  : 


fi  1*328 

S  =  Vll 


12 

For  very  small  M  and  low  R,  however,  Janour  has  shown  that  the  actual  drag 
coefficient  is  approximately  1,5  times  larger  than  that  given  by  Blasius. 
Hence  for  skin  friction  on  two  sides  of  the  dipole: 
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OkS 

for  M  <  0.2  (11  ) 

50  >  R  >  2 

where  R  is  based  on  the  dipole  dimension  vfhich  is  parallel  to  the  downward 
motion.  It  follows  from  equation  (2): 

M  <  0.2  (12) 

50  >  R  >  2 

If  the  dipole  falls  with  its  length  vertical  instead  of  width  vertical  then 
w  and  L  must  be  interchanged  in  the  above  equation, 

2.5  20  Km  altitude 

At  approximately  20  I<jn  altitude  it  is  found,  for  the  dipoles  discussed 

in  this  report,  that  the  Reynolds  number  falls  below  unity.  Such  a  low  number 
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is  outside  the  province  of  'aircraft*  aerodynamics  but  was  considered  by  Lamb 
in  his  theory  of  viscous  flow.  By  'viscous'  flow  is  meant  flow  in  which  a 
boundary  layer  of  air  molecules  is  permanently  attached  to  the  dipole  and  in 
which  the  mean  free  path  of  air  molecules  is  small  compared  with  any  dipole 
dimension.  The  flow  is  assumed  laminar  and  the  dipole  velocity  is  relatively 
low  (R  <  1  and  M/R  «  1 ) . 

Lamb's  theory  applies  to  a  cylinder  falling  with  its  length  horizontal 
and  unfortunately  is  obtained  in  a  form  which  contains  log  R  and  does  not 
permit  a  simple  solution  of  equation  (2)  for  the  velocity.  However  an  approxi¬ 
mation  to  Lamb's  equation,  correct  to  2C^6,  is^: 

S  <  R  <  2  .  (13) 

R 

The  above  formula  is  for  a  cylinder,  with  R  in  terms  of  the  cylinder  diameter 
(d),  Bairstow"'^,  however,  has  shovm  that  a  flat  strip  width  w  falling  with  its 
length  horizontal  has  the  same  motion  as  the  cylinder  provided 


d  =  0.83 w 

(strip  broadside- on  to  flow) 

(14) 

d  =  0,3 w 

( strip  edge-on  to  flow)  . 

(15) 

Substituting  in  (2)  it  is  then  found 


3 


7 


■  -1 . . 

U  j  *5  4  , 

\  /  5  |i  kp  w 


for  2  >  R  >  0.001  (16) 

M 

and  ~  <  0,01 


where  k  =  0,83  or  0,3  according  to  whether  to  strip  is  broadside-on  or  edge-on 
to  the  flow, 

2,4  Slip-flow  region 

As  the  altitude  is  further  increased  the  mean  free  path  of  air  molecules 

becomes  a  significant  fraction  (l-1C^)  of  the  boundary  layer  thickness  and  a 

phenomenon  known  as  *slip*  occurs;  the  layer  of  gas  adjacent  to  the  dipole 

begins  to  slide  at  a  finite  velocity  along  the  dipole  surface.  For  R  <  1 ,  the 
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effect  is  believed  to  occur  at: 

*  0.01  <  I  <  0.1  .  (17) 

The  only  available  formula  for  C„  in  this  region,  with  R  <  1 ,  appears  to  be 
1  c  ^ 

that  due  to  Tsein  for  a  horizontal  cylinder: 

(18) 


= 


8% 


R 


log^ 


1).  1,28  .  3.55 


This  is  a  modification  of  Lamb’s  formula*  for  finite  M/R.  Eqmtion  (l8)  is 
actually  Tsien's  formula  rearranged  (i)  so  that  R  is  based  on  the  cylinder 
diameter  instead  of  cylinder  radius  as  assumed  by  Tsien,  and  (ii)  for  the 
special  case  where  y  =  1 *41  and  where  air  molecules  are  assumed  to  reflect  in 
a  diffuse  way  on  striking  the  dipole.  The  formula  does  not  permit  an  explicit 
solution  for  V, 

9 

The  simplest  way  of  finding  V  appears  to  be  that  suggested  by  Barr  ;  by 
substituting  for  R  and  Cp  in  equation  (2)  it  is  found**: 


A.TtL.ii^.  ^  _  2  .S.8jl\ 

\  mg  •  0  J 

P  d(0.8  -  log^  r) 


(19) 


*  Lamb’s  equation  is  similar  to  equation  (l8)  for  M/R  =  0  except  that  Lamb 
has  (y~  2)  in  place  of  1.28. 

**  The  equation  given  in  Barr’s  report  omits  the  R  shown  in  the  numerator  of 

(19). 
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Thus  assuming  |i  and  o  are  constant,  it  is  possible  to  find  the  density  p  which 
corresponds  to  a  given  R,  Once  p  and  R  are  known,  it  follows  from  the  defini¬ 
tion  of  R  that: 


V  =  ^  ^  (20) 

d  p 

whilst  the  altitude  corresponding  to  V  may  be  deduced  from  altitude-density 
tables  (Table  1 )  . 

For  a  flat  strip,  it  is  assumed  equations  (I4)  and  (15)  still  apply. 

2.5  Free  molecule  flow  region 

At  greater  altitudes,  (above  70  Km)  the  mean  free  path  of  air  molecules 
becomes  much  greater  than  the  dipole  dimensions  (Table  2),  In  this  case  there 
is  no  boundary  layer  and  molecules  reflected  or  re-emitted  from  the  dipole 
surface  do  not  collide  with  air  molecules  until  they  have  travelled  a  consider¬ 
able  distance,  i.e,  the  dipole  does  not  cause  any  significant  distortion  of  the 
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air  stream.  This  situation  is  termed  ’free  molecule  flow'  and  drag  is 
determined  solely  by  molecule- surface  interaction,  being  independent  of  air 
viscosity.  It  has  been  shown  theoretically  that  free  molecule  flow  should 
occur  at 

M  >  3R  .  (21 ) 
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It  vdll  be  assumed  that  molecules  are  reflected  in  a  diffuse  manner  and 

that  the  dipole  temperature  is  equal  to  the  air  temperature.  It  is  then 
10 

found  that  for  a  cylinder  falling  with  its  axis  horizontal: 


0.5  -0.5s 

% _ e _ 


[(s  + 1.5) 


I^(0.5 


s^) 


+  (s  +0.5) 


1/0.5  s' 


)] 


% 


4s 


(22) 


where  I  and  I.  are  the  modified  Bessel  functions.  At  small  s,  I  ~  1  and 
0  1  o 

I^  ~  0,  The  quantity  s,  known  as  the  'molecular  speed  ratio’  is 

s  =  ~  (23) 


where  is  the  most  probable  molecular  velocity  ; 


C 


2 


ra 


2R'T 


W 


(24) 
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where  R’  is  the  gas  constant,  ¥  the  molecular  weight  of  air  and  T  the 
temperature. 

In  the  preliminary  calculations  it  was  noticed  that  for  window  dipoles 
below  100  Km, 

s  <  ^  (25) 

and  in  this  special  case  equation  (22)  simplified  to 


(26) 
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Barr  has  also  considered  the  case  of  the  cylinder  moving  with  its  axis  at 
right  angles  to  the  flow:  working  from  an  expression  for  which  is  more  general 
than  that  given  in  equation  (22)  he  found; 


0.5 


"D 


1  ,5  %  +  0.25  % 


i.5 


for  3  <  Omk 


(27) 


Evaluating  the  numerator,  this  is  found  to  be 


in  close  agreement  with  (26), 


Turning  to  the  flat  plate,  it  is  found 
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(28) 


0.5 

%  3 


e“(®  ®^“®)  +  ^0.5  g  0  A  +  erf(s  sin  6)  +  x  sin^0 


...  (29) 

This  equation  is  based  on  the  same  assumptions  as  (22).  0  is  the  angle  of 
attack,  i.e,  the  angle  between  the  surface  and  the  flow  direction  (vertical). 
erf(s  sin  0)  is  the  ’error  function’  and  at  small  s  sin  0; 


erf(s  sin  0) 


2 


% 


0.5 


s  sin  0 


(30) 


At  small  (s  sin  0),  equation  (29)  simplifies  to; 


10 


Oks 


c 


D 


~"q~V —  (1+4.142  sin^0) 
%  *  s 


Henoe  for  a  dipole  edge-on  and  length  either  horizontal  or  vertical 


Ldi 

s 


and  for  a  dipole  hroadside-on 


C 


D 


(31) 


(32) 


(33) 


Also  if  dipoles  are  in  random  orientations,  i.e.  the  probability  of  any  value 
of  0  is  proportional  to  sin  0,  then  the  mean  drag  coefficient  obtained  by 
integration  of  (3I )  over  all  0  is 


C 


D 


4.25 

s 


(34) 


which  is  nearly  the  same  as  the  value  for  the  horizontal  cylinder  (equation  (26)), 

As  in  the  previous  oases,  the  terminal  velocity  may  be  found  by  substitu¬ 
tion  of  Cjj  in  equation  (2).  For  example  (26)  yields: 


V 


lisN  -  1  ■ - 

j  *  2  p w 


(35) 


whilst  for  the  flat  dipole  broadside-on  to  the  air  flow  (53) 


V 


1 

2.9pw 


(36) 


3  LOW  ALTITUDE  MEASUREMENTS 

Although  a  number  of  raeas\arements  have  been  made  of  the  rate  of  fall  of 

window  dipoles  in  the  lower  atmosphere,  many  of  the  data  are  unfortunately 

rather  inaccurate  and  cannot  be  used  for  a  precise  check  of  the  theoretical 
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formulae.  Recently,  however,  Jiusto  and  Eadie  have  made  a  number  of  careful 
measurements  in  a  high-altitude  test-chamber,  simulating  altitudes  up  to  20  Km. 
Their  results,  together  with  details  of  the  dipoles,  are  reproduced  in  Pig.1. 
The  procedure  was  to  release  a  number  of  dipoles  and  then  note,  at  each 
simulated  altitude,  the  times  taken  by  the  fastest  and  slowest  dipoles  to  fall 
a  standard  distance  (2.5  metres).  A  record  was  also  kept  of  the  median  time. 
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Thus  the  three  curves  in  Fig.1  show  the  maximum,  median  and  minimum  velocity 
as  a  fimction  of  altitude# 

Also  plotted  on  Fig#!  are  points  calculated  with  the  aid  of  equations 
(l0),  (l2)  and  (16),  The  parameters  assumed  for  the  atmosphere  are  set  out  in 
Table  1,  It  can  be  seen  that  there  is  good  agreement  between  theory  and 
experiment;  the  slope  of  the  theoretical  curves  are  almost  identical  with  the 
measured  values  whilst  the  ratio  of  maximum  to  minimum  velocity  found  to  be 
6,8  from  the  experimental  data  compares  well  with  the  theoretical  values  of 
7.0  at  0  Km  and  7.4  at  20  Km, 

There  is  probably  little  significance  in  the  coincidence  of  the  median 
velocity  curve  with  the  theoretical  points  for  edge-on  horizontal  dipoles.  It 
is  commonly  observed  near  ground  level  that  most  dipoles  fall  with  their  length 
and  width  horizontal  so  presumably  the  median  res\ilted  from  a  few  dipoles  falling 
vertically  together  with  a  larger  number  falling  broadside-on.  Nevertheless  it 
is  useful  to  note  that  the  median  result  fits  the  edge-on  formula. 

The  fact  that  the  observed  maximum  velocities  are  slightly  less  than  the 
theoretical  values  might  be  explained  by  assuming  no  dipoles  fall  the  whole 
distance  with  their  lengths  exactly  vertical.  As  just  mentioned  the  most 
stable  dipole  orientation  appears  to  be  that  in  which  the  length  is  horizontal 
and  consequently  vertical  dipoles  might  tend  to  tiirn  and  so  have  reduced 
velocities. 

The  fact  that  no  dipoles  fall  quite  as  slowly  as  the  theoretical  minimum 
is  perhaps  rather  more  puzzling;  at  0  Km  the  difference  between  theory  and 
experiment  .is  approximately  15^,  The  authors  do  not  give  the  estimated  accuracy 
of  the  measurements  but  it  can  be  observed  that  a  dipole  falling  at  0,35  m/seo 
takes  7  seconds  to  fall  2,5  metres  and  hence  assximing  a  timing  to  an  accuracy 
of  0,2  second,  the  experimental  error  was  probably  no  greater  than  3?^,  i,e, 
only  one  fifth  of  the  discrepancy.  Possibly  no  dipole  fell  the  whole  distance 
with  its  surface  always  horizontal. 

The  authors  state  that  the  experimental  results  were  appropriately 
corrected  for  the  time  required  for  the  chaff  to  reach  terminal  velocity  but  no 
details  are  given  of  the  magnitude  of  the  correction. 

Measurements  were  also  made  with  S-band  dipoles.  These  had  the  same  width 
and  thickness  as  those  in  Fig.l  but  were  3*34  times  longer.  It  was  found  that 
the  rates  of  fall  for  the  longer  dipoles  were  the  same  as  those  in  Fig.l ;  in 
agreement  with  equations  (9),  (l0)  and  (l2)  which  all  state  that  V  is  propor¬ 
tional  only  to  the  mass  per  unit  length. 
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4  HIGH  ALTITUDE  MEAStlREKfEOTS 

The  most  oomprehensive  series  of  high  altitude  measurements  of  rate-of- 
fall  appear  to  be  those  reported  by  Smith  .  A  total  of  742  rate-of -fall  versus 
altitude  measurements  were  made  at  altitudes  in  the  range  40-84  Km  above  Nevada, 
the  Marshall  Islands,  Johnston  Island  and  V/hite  Sands,  In  each  case  the  rate- 
of-fall  was  measured  over  a  one-minute  interval.  The  curve  v;hich  best  fitted 
all  experimental  results  was  found  to  be: 

log^Q  f  =  5.014  -  2.G47  X  10“^  z  +  9.294  x  1o"^^  (37) 

where  f  =  rate  of  fall  is  ft/minute 
Z  =  altitude  in  ft. 

This  relation  is  shown  as  a  thick  curve  in  Pig, 2,  Also  in  Fig. 2  are  set  out 
details  of  the  dipoles  used  and  a  set  of  four  vertical  lines  which  indicate  the 
standard  deviation  in  individual  results  at  various  altitudes.  Smith  aotually 
gives  more  details  of  standard  deviations  but  only  four  representative  values 
are  included  in  Pig, 2  so  as  to  indicate  the  order  of  magnitude  without  obscuring 
other  detail. 

Some  of  the  Johnston  Island  measurements  discussed  by  Smith  were  also 
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analysed  in  greater  detail  by  Rapp  .  His  mean  curve,  for  altitudes  between 
72,5  and  80  Km  is  included  in  Fig. 2  and  can  be  seen  to  be  in  close  agreement 
with  equation  (37). 

As  the  dipoles  used  in  the  high  altitude  measurements  were  not  identical 
with  those  used  at  low  altitudes  (Section  3)>  the  results  in  Section  3  could 
not  be  included,  as  they  stood,  in  Fig. 2,  In  view,  however,  of  the  good  agree¬ 
ment  between  theory  and  experiment  at  low  altitudes  it  was  felt  useful  to  scale- 
up  the  median  low-altitude  results,  with  the  aid  of  equation  (12)  and  plot  these 
also  on  Fig. 2,  The  resulting  curve  is  shown  as  a  thick  line  between  0  and 
20  Km,  An  interpolation  has  been  made  in  Pig. 2  between  the  low  altitude  results 
(below  20  Km)  and  the  high  altitude  results  (above  40  Km)  in  order  to  gain  some 
ideas  of  the  likely  variation  in  V  at  all  altitudes  below  83  Km. 

Also  drawn  on  Fig. 2  are  the  theoretical  curves  for  the  slowest  falling 
dipoles  (w  and  L  in  a  horizontal  plane).  Each  curve  is  drawn  as  a  solid  line 
for  those  altitudes  at  which  R  and  M/R  are  within  the  theoretical  boundary 
conditions  set  out  in  Section  2  and  as  a  dashed  line  at  other  altitudes.  The 
'free-moleoule  theory’  curve  is  not  dravm  above  90  Km  altitude  because  the 
molecular  speed  ratio  s  then  exceeds  the  value  0,5  assumed  in  the  approximate 
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equations  (32)-(36),  In  the  free-moleoule  zone,  a  second  theoretical  curve 
(dotted  and  labelled  *LtAX*)  has  been  drawn-in  for  the  fastest  falling  dipoles 
(equation  (32)). 

Comparing  theory  and  experiment  it  may  be  observed  that  at  each  altitude 
there  is  remarkably  good  agreement  between  the  slopes  of  the  experimental  curve 
and  the  slope  of  the  appropriate  theoretical  curve.  The  four  theoretical  curves 
form  a  continuous  S-shape  very  similar  to  the  experimental  result. 

In  Fig. 2  it  can  be  seen  that  the  experimental  curve  takes  a  fairly  sharp 

bend  upwards,  i.e,  the  slope  departs  from  that  predicted  by  the  viscous  drag 
formula,  at  an  altitude  of  about  65  Km.  This  is  almost  exactly  the  altitude 
at  which  the  mean  free  path  of  air  molecules  becomes  equal  to  the  dipole  width 
(Table  2),  Here  the  ratio  M/R  ~  0.7.  At  higher  altitudes  the  slip-flow 
formula  appears  to  be  in  reasonable  agreement  with  experiment  at  least  until 
M/R  >  3(75-80  Km)  where  the  free-molecule  formtila  applies,  (incidentally  the 
slip  flow  formula  (19)  cannot  be  applied  at  very  low  Reynolds  numbers  because 
the  numerator  becomes  negative.) 

One  unexpected  result  is  that  the  viscous  drag  and  slip-flow  formulae 
both  seem  to  apply  at  higher  altitudes  than  predicted.  According  to  theory^^, 
slip  flow  should  commence  and  the  viscous  drag  formula  should  begin  to  fail  at 
M/R  ~  0,01,  (i.e,  35  Km)  o.f,  the  experimental  value,  M/R  ~  0,7  (65  Km).  Further, 
whereas  the  slip-flow  regime  is  normally  considered  to  end  at  M/R  ~  0,1  (50  Kin), 
in  practice  the  equation  appears  reasonably  accurate  up  to  U/R  ~  3  (75  Km),  In 
other  wcrds  significant  changes  in  the  drag  coefficient,  due  to  finite  molecular 
mean  free  path,  appear  to  occur  only  at  much  larger  mean  free  paths  than 
estimated  theoretically, 

10 

A  further  point  of  interest  is  that  in  a  discussion  of  the  theory,  Emmons 
(1958)  states  (i)  that  no  experimental  results  are  available  to  check  the  slip- 
flow  formula  for  the  drag  on  a  cylinder  (eqmtion  (l8))  and  (ii),  in  the  region 

0.1  <  f  <  3 

called  the  'transition  regime*,  veiy  little  is  known  about  either 
theoretically  or  from  experiment*.  The  present  results  suggest  that  equation  (l8) 

*  Some  laboratory  measurements  have  been  made'  in  the  transition  regime  for 
large  Mach  numbers  (M  >  2)  but  apparently  none  have  been  made  for  the  low  Mach 
numbers  encountered  in  the  present  study. 
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is  in  fact  correct  but  that  the  slip-flow  region  extends  right  up  to  the  free- 
molecule  zone,  i.e.  there  is  no  need  to  define  a  separate  Hransition  regime’. 

The  theoretical  curves  in  Fig, 2  are  for  the  slowest  falling  dipoles. 
Presumably  the  median  experimental  results  are  approximately  1,3  times  greater 
because  many  dipoles  tend  to  fall  edge-on  rather  than  broadside-on.  At  the 
higher  altitudes  the  curves  for  minimum  velocity  and  for  maximum  velocity  run 
at  approximately  1,5  standard  deviations  either  side  of  the  mean  experimental 
curve:  talcing  into  account  the  fact  that  at  the  highest  altitude  the  error  in 
a  single  measurement  of  rate-of-fall  may  be  of  the  order  of  10  m/sec  (Rapp'), 
it  is  probable  that  most  of  the  scatter  in  the  experimental  data  arose  from 
measurement  errors  and  there  is,  therefore,  no  strong  evidence  that  any  dipoles 
fell  with  velocities  outside  the  theoretical  limits. 

One  assumption  made  in  the  foregoing  discussion  is  that  the  measurements 

were  all  made  on  dipoles  falling  at  the  terminal  velocity.  This  point  was 

7 

discussed  by  Rapp  and  was  also  considered,  for  the  free  molecule  region,  by 
9 

Barr  .  Both  conclude  that  the  assmption  is  very  probably  correct. 

At  low  altitudes  it  was  found  that  the  median  results  agreed  closely  with 
the  theoretical  values  for  horizontal,  edge-on  dipoles.  To  see  if  this  agree¬ 
ment  also  held  at  greater  altitudes,  the  experimental  results  in  Fig. 2  were 
reproduced  in  Pig, 3  together  with  the  theoretical  curves  for  dipoles  falling 
with  the  width  vertical,  (The  free-molecule  curve  was  obtained  by  fitting  to 
the  slip-flow  curve  at  the  point  where  M/R  =  3  corresponds  to  Cp  ~  2,6/s,) 

Fig, 3  shows  that  the  ’edge-on'  theoretical  curves  are  in  close  agreement 
with  experiment  and  it  appears  that  the  relation:  (median  velocity)  ~ 
(theoretical  velocity  for  horizontal,  edge-on  dipoles)  holds  at  all  allitudes 
investigated, 

5  DISCUSSION 

5*'1  Equation  for  rate-of-fall 

It  will  have  been  noticed  that  the  theory  outlined  in  Section  2  is 
complicated  and  the  individual  equations,  particularly  those  for  slip-flow,  are 
tedious  to  evaluate.  No  less  than  four  separate  theories  are  employed  and 
unfortunately,  vintil  numerical  values  of  V  and  R  and  M  have  been  worked-out  for 
a  given  altitude  it  is  not  possible  to  decide  whether  the  particular  equation 
used  was  in  fact  the  correct  one.  It  would  clearly  be  useful  to  have  an 
approximate  formula  valid  over  the  whole  range  of  altitude. 
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To  obtain  such  an  approximate  equation  it  may  be  noticed  that  the 
viscous  drag  theory  is  only  slightly  in  error  at  all  altitudes  up  to  that  at 
which  there  are  pronounced  slip  effects.  Near  ground  level  the  error  in  the 
viscous  drag  theory  is  small  and  would  be  less  with  thinner  dipoles  (smaller  R). 
Also,  onoe  slip-flow  effects  have  become  important,  only  a  small  increase  in 
altitude  is  required  before  entering  the  free  molecule  theory  domain.  Thus  an 
approximate  theory  might  be  based  solely  on  viscous  drag  and  free  molecule 
theories. 

One  way  of  combining  the  two  theories  is  simply  to  add  the  two  velooities 
together; 


Here  the  first  term  on  the  right  hand  side  is  the  viscous  drag  term  (equation 
(16))  with  |i  set  at  the  mean  value  1,5  x  10”^  ogs  units,  whilst  the  second  term 
is  the  free  molecule  equation  (3I )  with  =  4  x  10^  and  =  2,6/s  (Section  4), 
At  low  altitudes  the  second  term  is  negligible  and  conversely  at  high  altitudes 
the  first  term  oan  be  ignored. 

Substituting  for  g  in  equation  (58)  and  re-arranging  into  a  more 
convenient  form: 


35000  m  0.02  m 
_ o_  ^  _ o 

6  Vm^  p  w  P  w 


vfhere  V 
m 


P 


terminal  velocity  -  cms/second 

mass  per  unit  length  (m^  =  m/b)  gms/om 

dipole  width  (cm) 

atmospheric  density  in  ogs  units  (Table  1 ) . 


(39) 


This  equation  is  the  approximation  for  V  for  all  altitudes  below  100  Km, 

Points  calculated  from  (39)  are  shown  as  small  crosses  in  Pig. 3  and  can  be  seen 
to  lie  close  to  the  experimental  curve,  the  error  being  less  than  the  standard 
deviations  in  the  measurements  (Pig,2), 


5,2  Slowest-falling  dipoles 


In  measurements  of  wind  intensity  and  direction,  it  is  desirable  to  use 
slow-falling  dipoles  for  two  separate  reasons. 
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(i)  to  reveal  detail  in  the  wind  structure,  e.g.  measurements  with  fast- 
falling  copper  dipoles  revealed  considerably  less  detail  than  those  with  slower- 
falling  nylon  dipoles^, 

(ii)  to  extend  the  useful  lifetime  of  each  window  cloud. 

At  the  lower  altitudes,  it  is  clear  from  equation  (39)  or  from  the  theory 
in  Section  2,3  that  apart  from  very  slight  variations  due  to  the  sixth  root 
term,  the  terminal  velocity  varies  only  with  the  mass  per  unit  length; 


V  oc  m 

0 


Some  types  of  chaff  mentioned  in  the  literature 


5,6,9 


are  as  follows; 


Type 

ra 

gms 

L 

oms 

m 

0 

aluminium  strip  (Fig.1 ) 

X  10"^ 

0.95 

1 .52 

X  10"^ 

0.62 

aluminium  strip  (Fig.2)* 

7.5 

5.0 

1.5 

nylon  cylinder  (0.0035") 

8.1 

4.8 

1.7 

nylon  cylinder  (0.008") 

19.5 

A«8 

4.1 

nylon  cylinder  (0,012")*'" 

47 .6 

4.8 

10 

copper  cylinder  (O.OlO")* 

216 

4.8 

45 

(*  denotes  types  reported  used  in  atmospheric  wind  measurements) 

It  can  be  seen  that  the  0.012  inch  nylon  used  in  wind  measurements  was 
relatively  fast-filling.  The  strips  or  thinner  nylon  wo’xLd  have  been  preferable. 

At  the  highest  altitudes  the  drag  coefficient  is  approximately  the  same 
for  strips  and  cylindrical  dipoles  (equations  (34)  and  (26))  and 

m  m 

V  oc  — ^  or  .  (A1  ) 

Thus  for  a  strip  where  m^  =  .  wt; 

V  oc  .  t  {k2) 

where  p^  is  the  density  of  the  dipole  material. 
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Similarly  for  a  cylindrical  dipole  where  =  Pjq  • 

V  oc  ,  7i:d/4  ,  (43) 

It  follows  that  for  minimum  rate-of-fall  at  high  altitude,  the  material  density 
should  be  kept  low  and  also  strips  are  probably  preferable  to  cylinders  because 
the  thickness  (t)  of  a  strip  can  probably  be  made  less  than  the  diameter  (d)  of 
a  metal  or  nylon  thread.  Typical  values  are  t  =  1,14  x  lO"'^  cms  (strip  dipoles 
in  Fig.l )  and  d  =  8,89  x  10”^  cms,  (thinnest  nylon  dipoles  discussed  by  Barr^), 

6  CONCLUSIONS 

(1)  Measurements  of  the  rate~of-fall  of  window  dipoles  at  high  altitudes 
appear  to  be  in  good  agreement  with  the  theory.  At  low  altitudes  viscous  drag 
effects  are  dominant  whilst  above  75  free-raoleoule  theory  applies.  In  the 
intermediate  'slip-flow*  region,  the  experimental  results  appear  to  confirm 
Tsien's  equation  for  the  drag  on  a  cylinder. 

(2)  At  low  altitudes  (e.g.  below  50  Km)  the  slowest  falling  dipoles  are  those 
which  have  lowest  mass  per  unit  length.  At  high  altitudes,  (e.g.  above  70  I^m) 
the  terminal  velocity  is  proportional  to  the  density  and  thickness  of  the 
metallic  strip, 

(3)  In  view  of  the  possible  number  of  ways  of  falling  the  agreement  between 
the  theory  and  practice  examined  in  this  report  might  seem  to  be  fictitious; 
there  is  no  evidence  available  to  show  what  would  happen  for  dipoles  of 
different  geometries  or  at  radically  different  altitudes. 


Table  1 

ATMOSPHERIC  PRESSURE  /  DENSITY.  VISCOSITY  AM)  TEMPERATURE 


All  quantities  are  expressed  in  ogs  units 


Altitude 

Pressure  P 

Density  p 

Viscosity  p 

Temperature  T 

Km 

0 

1  X  10^ 

1,2  X  10"^ 

1.8  X  10"^ 

°K 

288 

10 

2.6  X  10^ 

4x10^ 

1 .4  X  10“^ 

223 

20 

5.4  X  10^ 

9  X  10"^ 

1 .4  X  10“^'' 

217 

30 

1.2  X  10^ 

2  X  10~^ 

1 .5  X  10"^ 

227 

AO 

2.9  X  10^ 

4  X  10"^ 

1 .7  X  10~^ 

250 

30 

800 

1  X  10"^ 

1 .7  X  10~^ 

270 

6o 

210 

3  X  10“^ 

1.7  X  10“^ 

250 

70 

54 

9  X  10“® 

1.4  X  10“^ 

220 

80 

10 

2  X  10“® 

1 .2  X  10"^ 

180 

90 

1 .6 

3.2  X  10“^ 

1.2  X  10"^ 

180 

100 

0.29 

5  X  lO"'"^ 

! 

210 

19 


Table  2 


APPROXIMATE  REYNOLDS  MJMBERS.  MACH  NUMBERS  AM)  MEM  FREE  PATHS 

The  numbers  are  oaloulated  from  the  data  in  Table  1 
and  the  observed  median  velocities  (Fig,2).  The 
Reynolds  nmbers  are  based  on  the  dipole  width. 


I 

Altitude 

R 

M 

m/r 

mean  free  path 

■ 

20 

0.002 

0.0001 

cms 

7  X  10"^ 

15 

0,004 

0.0003 

2  X  10”^ 

5 

0.007 

0.0014 

9  X  10”^ 

1.5 

0.009 

0.006 

4.4  X  10"^ 

40 

0.3 

0.01 

0.03 

0,002 

50 

0.1 

0.012 

0.1 

0.008 

60 

0,04 

0.017 

0.5 

0.026 

70 

0.025 

0.03 

1 .4 

0.093 

80 

0.025 

0.13 

5 

0.41 

90 

0.02 

0.7 

30 

2.5 

100 

0.02 

3.5 

200 

16 

Ti;onfnC.':L 

FLi)'; 


HB-EIET 
HI  3 


adeed: 


GP.OUHB,  MB. 


SYMBOLS 


A  dipole  projected  area  normal  to  airflow  (sq  cms) 

0  velocity  of  sound  (equation  (5)  -  approximately  3x10^  oms/seo) 
Cp  drag  coefficient 

C  most  probable  molecular  velocity  (equation  (24)  -  approximately 
(4  X  10^  cms/sec) 

d  diameter  of  cylindrical  dipole  (cms) 
f  rate  of  fall  in  vmits  of  feet  per  minute 

2 

g  acceleration  due  to  gravity  (98I  cms/sec  at  ground  level,  957  cm 
80  Km) 

L  dipole  length  (cms) 

m  dipole  mass  (gms) 

m^  dipole  mass  per  unit  (gras/cm) 

M  Mach  number  (equation  (4)) 

P  air  pressure 

R  Reynolds  number  (equation  (3)) 

R*  gas  constant 

s  moleoular  speed  ratio  (equation  (23)) 
t  dipole  thickness  (cms) 

T  absolute  temperature  (Table  1 ) 

V  terminal  fall  velocity  through  still  air  (cms  per  sec) 
w  dipole  width  (cms) 

W  molecular  weight  of  air 
Z  altitude  in  feet 

Y  ratio  of  specific  heats  of  air  molecules  (1 .41 ) 
p  air  density  (Table  1 ) 

p^  mean  density  of  dipole  material 
p.  air  viscosity  (Table  1 ) 
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FIG.  I.  VELOCITIES  AT  LOW  ALTITUDES 
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FIG.  3.  VELOCITIES  AT  HIGH  ALTITUDES  - 
BEST  THEORETICAL  ESTIMATES 
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